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Abstract 

Four-particle correlations in fermion systems at finite temperatures are in- 
vestigated with special attention to the formation of a condensate. Instead 
of the instability of the normal state with respect to the onset of pairing de- 
scribed by the Gorkov equation, a new equation is obtained which describes 
the onset of quartetting. Within a model calculation for symmetric nuclear 
matter, we find that below a critical density, the four-particle condensation 
(a- like quartetting) is favored over deuteron condensation (triplet pairing). 
This pairing-quartetting competition is expected to be a general feature of 
interacting fermion systems, such as the excition-biexciton system in excited 
semiconductors. Possible experimental consequences are pointed out. 
Keywords: Bose-Einstein condensation, superfluidity, strongly coupled sys- 
tems, nuclear matter, a-particle matter, exciton-biexciton system. 
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One of the most amazing phenomena in quantum many-particle systems is the formation 
of quantum condensates. Of particular interest are strongly coupled fermion systems where 
bound states arise. In the low-density limit, where even-number fermionic bound states 
can be considered as bosons, Bose-Einstein condensation may be expected to occur at low 
temperatures. At present, condensates are investigated in systems where the cross-over 
from Bardeen-Cooper-Schrieffer (BCS) pairing to Bose-Einstein condensation (BEC) can 
be observed, see Q. Among very different quantum systems such as electron- hole exciton 
system in semiconductors, atoms in traps at extremely low temperatures, etc., nuclear matter 
is particularly suited to study correlation effects in a quantum liquid. 

An indication of strong correlations in nuclear matter is the formation of bound states. 
The interaction in the singlet (S = 0) channel is not strong enough to form a bound state, 
whereas in the neutron (n) - proton (p) triplet (S = 1) channel a bound state, the deuteron, 
arises. In nuclear matter, with increasing density the bound states are modified and will 
disappear at the so-called Mott density 0. In particular, treating the two-particle Green 
function in ladder Hartree-Fock approximation, an effective wave equation (in matrix nota- 
tion) ipx = K 2 (Ex) ipx for the quantum state A can be derived. Explicitly this reads 

Va(12) = ^K 2 (12,1'2',^a)^a(1'2') (1) 

1'2' 

with 

KM2, 1'2', z) = V(12, 1'2') 1 - - . (2) 
2V ' ' ; V ' J z-E(l)-E(2) K ' 

The influence of the medium is contained in the single-particle energy E(l) = p\/2m + 
£ 2 1/(12, 12) cx /(2) and in the Pauli blocking term [1 - /(l) - f(2)]. Here /(l) = 
[exp{E(l)/T — fi/T} + is the Fermi distribution function and T' stands for momentum, 
spin, and isospin coordinates, whereas V(12, 1'2') is the antisymmetrized matrix element of 
the two-body interaction. 

Including medium modifications of two-particle states, a generalized Beth-Uhlenbeck for- 
mula N/Qq = ni(T, //) + n corr (T,fi) has been derived UK], where the uncorrelated density 
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rii(T, //) = Qq J2i /(I) is given by the quasi-particle contribution (fl being the normaliza- 
tion volume). The correlated density n corr (T, //) is obtained from the two-particle T matrix. 
This approach has been widely applied to ionic plasmas as well as to the electron-hole ex- 
citon system in excited semiconductors ||. Taking into account the Mott effect, general 
features of the composition of nuclear matter as function of density and temperature are 
given in |J. 

At low temperatures it is well known that in nuclei, nuclear matter and neutron matter 
(neutron stars) superffuidy can arise in the singlet channel [|7|-f| • A theoretical description of 
this superfluidity can be achieved by treating the Gorkov equation ip 2 = K 2 (^1+^2) 1P2 for the 
critical temperature T 2 C as a function of the chemical potential. It allows the determination 
of T s c or T 4 C for the singlet and triplet channels, respectively. 

The solution of the Gorkov equation has been considered by different authors using 
realistic nucleon-nucleon interactions. It has been found that in comparison with the singlet 
channel, in the triplet channel the transition to superfluidity should arise at relatively high 
temperatures ||T0| , |TTf . This is a consequence of the stronger interaction in the triplet channel 
which leads to the formation of a bound state (deuteron) in the low-density limit where 
/ -C 1. Estimates give a value of the critical temperature up to T t c m 5 MeV at one third of 
the nuclear matter density. At the same time, at zero temperature a large gap arises [ H \ . An 
interesting feature of the triplet pairing in symmetric nuclear matter is the cross-over from 
Bose-Einstein condensation of deuterons at low densities to BCS neutron-proton pairing at 
high densities ||. 

In spite of the relatively strong interaction, cf. also calculations with effective pairing 
forces by Goodman |12|] , triplet pairing seems less apparent in nuclear structure systematics 
T3fl . However, it should become important for heavier N = Z nuclei produced in the new 



radioactive beam facilities. 

In this letter we show that in the low-density region the transition to triplet pairing 
is not realized, because four-particle correlations are more dominant there. Obviously, at 
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chemical equilibrium, in the low-density region at low temperatures the dominant part 
of nuclear matter will be found in a-particles which are much stronger bound than the 
deuteron. Therefore, the triplet pairing (Bose condensation of deuterons) has to compete 
with quartetting (Bose condensation of a-particles). 

The four-particle correlations are obtained from the four-particle Green function (com- 
pare ||14|| ) that is given in ladder Hartree-Fock approximation by the equation 

G 4 (1234, 1W4', z) = /(1)/ ' 2)/(3 » /(4) *'" fe t'^ 
4V ' ; £ 4 (1234) z-E 4 (1234) 

+ Y, ^4(1234, 1"2"3"4", ^)G 4 (1"2"3"4", l'2'3'4', z) , (3) 

1"2"3"4" 

where we use the abbreviation E n (12 . . . n) = E(l) + E(2) + ■ • • + E(n), and g n (12 . . . n) = 
[exp(£' n (12 . . . n) — nfi)/T — being the Bose distribution function. 

The instantaneous part of interaction kernel is obtained by using the technique of Mat- 
subara Green functions as 

K 4 (1234, 1'2'3'4', z) = V(U, l'2* /(1)/(2) 



g 2 {\2) z-E 4 (1234) 
+ perm. , (4) 

where the terms obtained by renumbering are not given explicitly. We have used the identity 
/(l)/(2) • • • f(n) - /(l)/(2) • • • f(n) = g~\\2 . . . n)/(l)/(2) • • • f(n) with / = 1 - /. 

Near a pole the Green function G 4 can be factorized, G 4 (1234, r'2'3'4', z) ~ 
^(1234)^*(l'2'3'4')/(-2 - E v ). The eigenvalues E v and eigenstates ^(1234) follow from the 
solution of the four-particle Schrodinger-like wave equation ip u = K^(E V ) ip u . The eigenvalue 
Eo(T,fi) of the lowest bound state (a-particle) depends on temperature T and chemical po- 
tential [L and also on its center-of-mass momentum due to the medium- dependent self-energy 
as well as the phase space occupation factors. 

From the four-particle Green function, the four-particle density matrix is obtained as 

< CLi CL2 a t a t Q4'03'G2'Gl' > 

= / — 2 4 (cu)ImG 4 (1234, r'2'3'4', u - iO) . (5) 

J 7T 



Obviously, the four-particle density (||) diverges when, at a given temperature T, the chem- 
ical potential takes the value \i c = E (T, /x c )/4. Then, the delta function produced by the 
pole of G4 coincides with the singularity of the Bose function. Values for the chemical po- 
tential fi exceeding the value of the lowest bound state energy E (T,fi) are not admissible 
because the diagonal elements of the four-particle density are positive definite. 

The main objective of this letter is to give an estimate of the critical temperature T£(fi) 
for the onset of a four-particle condensate. This is obtained by solving the equation 

^ 4 (1234)= ^4(1234, r'2'3'4', 4/i)^ 4 (l'2'3'4') . (6) 

1'2'3'4' 

In the low-density limit, where the distribution functions occurring in K4 are small compared 
with 1, this equation coincides with the Schrodinger equation for the four-particle bound 
state, i.e. the a-particle in free space. 

To discuss the competition between two-particle pairing and the condensation of 
four-particle states, we perform an exploratory calculation for a simple model system 
which contains the formation of a two-particle bound state in the triplet channel and 
a singlet four-particle bound state. We use a separable interaction of Gaussian form 
V s /t(12, 1/2') = — \s/t^o 1 exp(— P12/& 2 ) exp(— p'^/fo 2 ) <5(gi2 — q[ 2 ) with the relative momen- 
tum pu = (pi —p2)j1 and the center-of-mass momentum qyi = Pi + P2- At given range 
parameter b = 1.54 fm~ , the interaction strengths are adopted as \ t = 1213.8 Me Vfm 3 and 
A s = 536 Me Vfm 3 to reproduce the free deuteron binding energy E t = — 2.225 MeV in the 
neutron-proton triplet channel as well as the free a binding energy E = —28.29 MeV for the 
four-particle ground state (see below). In principle it is possible to extend the calculations 
to more realistic nucleon-nucleon potentials as given in ||15|| . 

We begin with the discussion of neutron-proton triplet pairing. The critical temperature 
T t c obtained from the solution of the Gorkov equation is presented as a function of the reduced 
chemical potential /z* in Fig. [I] and as a function of the uncorrelated density rii (T, /i) = 
^0 1 Si /(l) in Fig. (dotted lines). For the sake of simplicity the self-energy shift is 
taken to be a constant that is incorporated into a shift of the chemical potential, /i* = 
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/i-E 2 ^(12,12)ex/(2), at Pl = 0. 

The four-particle wave equation (|3|) is solved within a variational ansatz where the four- 
particle wave function with zero total momentum is given as a product of two wave functions 
for the relative motion (quasi deuterons) <f)(p) and the center-of-mass wave function ip(q), 

^ 4 (1234) = 0(^12)0(^34) ^(912)^12+534,0 • (7) 

The wave function for the relative motion 0(p) is optimized in a space of functions containing 
the exact solution for the wave function of the in- medium deuteron (Eqs. (|l],[|)), whereas for 
the center-of-mass wave function ifi(q) a simple Gaussian form is adopted. 

The results for the critical temperature T 4 as a function of the chemical potential \i* 
and the uncorrelated density n are displayed in Figs, p] and Q, respectively (solid lines). In 
the low-temperature limit, with increasing chemical potential the transition to quartetting 
occurs prior to the pairing transition. This is a consequence of the fact that the value Eq/4 
for the four-particle bound state lies below the value E t /2 of the triplet bound state. The 
Gorkov equation, which predicts a normal state for temperatures T > T t c , is not applicable 
in that region, because already at a higher temperature T| the normal state was removed 
due to the onset of quartetting. 

A simple argument for the behavior of T|, T t c as a function of n\ in the limit of 
low densities can be given from the law of mass action. Neglecting medium correc- 
tions the critical temperatures T t c and T 4 as functions of n\ are found as solutions of 
Hi = 4:(mT t c /2irh 2 ) 3 / 2 exp{E t /2T t c } and n x = 4(mT 4 c /27rft 2 ) 3 / 2 exp{£ /4T 4 c }, respectively. 

With increasing density the critical temperature T 4 approaches T t c from above at a critical 
density n\ = 0.03/m™ 3 (/iQ ps lOMeV). For ri\ < n\ (/x* < /ig) the four-particle condensate 
arises whereas for rii > n\ the system goes over to a BCS pairing state. 

The fact that at higher densities a-particle condensation disappears faster than deuteron 
Cooper pairing is at first sight quite astonishing because of the very strong binding of a's 
in free space. This feature, one of the main results of our work, can, however, be explained 
qualitatively with simple arguments and is generic for any Fermi system where pairing and 
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quartetting may interfere. Indeed the competition between pair and quartet condensation in 
attractive Fermi liquids is known also from the field of semiconductors, where both exciton 
and biexciton condensates have been suggested. The qualitative features of that competition 
are easily understood. Since they do not depend crucially on the triplet nature of the 
condensate we discuss the simpler case of singlet pairing as a function of density. In the 
weak coupling limit (high density), the ground state is described by the standard BCS wave 
function. However weak the attraction is, a pairing amplitude (aa) appears due to the finite 
density of states at Fermi level for pairs with zero total momentum q (the two-particle kernel 
has a logarithmic singularity). In contrast a quartet condensate (aaaa) ^ does not exist 
on its own. Due to the exclusion principle the corresponding density of states for q = 
vanishes at Fermi level and it takes a minimum attraction to develop a Cooper pole. 

In the opposite limit of dilute systems with attraction strong enough to bind individual 
pairs and quartets, an atomic regime prevails. Single pairs (aa) may be viewed as bosons ($) 
and the quartet as a bound pair of two bosons. The competition between particle and pair 



condensation in Bose liquids was discussed in [T6j . In the dilute limit we expect "molecular" 



Bose condensation of bound bosons, as described long ago by Valatin and Butler [|17|. In 



the dense limit the quartets dissociate and single boson condensation takes over. As the 
density decreases, a variational Ansatz [l6j shows that the condensate fraction uq = \(&)\ 2 



goes down, vanishing at a critical n c beyond which only ($$) survives (n c also marks the 
appearance of a gap). 

Altogether, we expect a quartet atomic condensate for very low densities and a BCS pair 
condensate for very high densities. In between the interpolation depends on the pair and 
quartet binding energies, E 2 and 2E 2 + E±, as compared to the Fermi energy. If E 2 ^> E4 
the pair condensate enters the atomic regime well before it disappears in favor of quartets. 
If on the other hand E 2 < E 4 the system goes directly from a BCS pair superfluid to a 
quartet Bose-Einstein condensation. 

Possible signatures of quartetting have been discussed in the past in the context of nuclear 



structure [flql . According to the present studies quartet correlations may be of importance 



in the outer regions of nuclei. We can make a rough estimate for the critical temperature 
of a - particle condensation in nuclei in converting T£(n) of Fig. 2 via the Local Density 
Approximation into a radius dependence T£(r) and averaging over a typical nuclear density. 
This procedure applied to ordinary neutron-neutron pairing yields T c w 1 MeV in good 
agreement with more microscopic calculations [113. In this way we predict that the increase 
of the critical temperature due to a - particle condensation is about 0.2 MeV. In spite of our 
limited variational approach, we think that this value should be correct within 30 to 40%. 
Also the critical density n\ = 0.03 fm~ 3 for the onset of quartetting should be reliable within 
these limits and therefore our estimate can be of interest for clustering phenomena in stellar 
collapse and the outer crust of neutron stars 0. We believe that in the far tail of nuclear 
densities (n < n^) quartetting can occur as well as in a number of nuclear processes such 



as the emission of a-particles from the heaviest elements |2(J an d unusually large a-decay 
rates of neutron deficient lead isotopes jH]]. Bose- Einstein effects may also be of relevance 
in heavy-ion collisions at a relatively late stage where the temperature of expanding matter 
has dropped below 5 MeV. 

In conclusion, we have shown that prior to the transition to pairing, a transition from 
the normal state to four-particle condensation according to Eq. (^j can occur. However, 
as a generic phenomenon 4-body condensation always looses against weak coupling Cooper 
pairing at sufficiently high densities. For demonstration purposes we have considered nu- 
clear matter as a strongly coupled quantum liquid. The possibility of quartetting and higher 
correlated condensates is of interest also in other systems, e.g., electron- hole pairs in semi- 
conductors forming bi-excitons [f[|,p2] or quarks forming two-pion states |23] . 

G. R. and A. S. appreciate useful discussions with G. Bertsch. P. S. is grateful to 
Deutsche Forschungsgemeinschaft (DFG) for financial support and to University of Rostock 
for pleasant hospitality. 
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FIG. 1. Critical temperatures for the onset of quantum condensation in symmetric nuclear 
matter, model calculation. The critical temperature of the onset of two-particle pairing T t c is 
compared with T| for the onset of a four-particle condensate, as a function of the chemical potential. 



10 



4r 
10 



T t ° (deuteron pairing) 
T 4 ° (a-quartetting) 




10 10 10" 

density n, [frrf 3 ] 



10 u 



FIG. 2. The same as Fig. but as a function of the uncorrelated density n\ . 
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